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Abstract: The Gallavotti-Cohen fluctuation theorem suggests a general symmetry in 
the fluctuations of the entropy production, a basic concept in the theory of irreversible 
processes, based on results in the theory of strongly chaotic maps. We study this 
symmetry for some standard models of nonequilibrium steady states. We give a general 
strategy to derive a local fluctuation theorem exploiting the Gibbsian features of the 
stationary space-time distribution. This is applied to spin flip processes and to the 
asymmetric exclusion process. 

Dedicated in honor of Robert Minlos on the occasion of his 70th birthday. 



1 Introduction 

A basic feature of equilibrium systems is that the restriction to a subsystem is again 
in equilibrium and with respect to the same microscopic interaction, for the same 
temperature, pressure and chemical potential. We can imagine cutting out a much 
smaller but still macroscopic region from our system and we will still find the same 
equilibrium state apart from possible boundary effects. Mathematically, this is ex- 
pressed via the DLR-equation stating that the local conditional probabilities of a 
Gibbs measure coincide with the corresponding finite volume Gibbs measures. This 
really amounts to the fact that, for Gibbs measures, the ratio of probabilities for 
two different microscopic configurations that are identical outside a finite volume, is 
given by the Boltzmann factor ~ exp (3 AH for relative Hamiltonian AH depending 
continuously on the configuration far out. In other words, relative energies make 
sense and they can be written as a sum of sufficiently local interaction potentials. 
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Important consequences are found in the theory of equihbrium fluctuations and in 
the framework of the theory of large deviations; Robert Minlos was among the very 
first to develop a mathematical theory of this Gibbs formalism. In particular in 
the description of the thermodynamic limit of Gibbs measures is given. We indeed 
usually have in mind here (very large) spatially extended systems which are moni- 
tored locally. 

The idea that time does not enter this equilibrium description is further stimulated 
from the fact that at least in classical statistical mechanics, the momenta (entering 
only in the kinetic energy part of the Hamiltonian) can be integrated out at once from 
the partition function. Time enters already more explicitly in equilibrium dissipative 
dynamics such as via Langevin equations (Markov diffusion processes) or Glauber 
dynamics (for spin relaxation). Yet, under the condition of detailed balance, the sta- 
tionary dynamics is microscopically reversible and the past cannot be distinguished 
from the future. The equilibrium steady state probability distribution on the space- 
time histories still has a Gibbsian structure with as extra bonus that the restriction to 
a spatial layer (at fixed time) is still explicitly Gibbsian and in fact, the restriction of 
the dynamics to a subregion still satisfies the detailed balance condition with respect 
to it. It was again Robert Minlos who was among the pioneers in the subject of 
space-time Gibbs measures and who, with his great experience in cluster expansion 
techniques and his love for field theory, saw the advantage of the space-time approach 
in the construction of solutions to infinite dimensional Markov diffusion processes, 
see e.g. ||. 

While the condition of detailed balance reflects a symmetry (the space-time distribu- 
tion is time- reversal invariant), the fact that the space-time distribution is Gibbsian 
(at least in some sense) does however not at all depend on it. In other words, the 
fact itself that the spatio-temporal probability distribution enjoys Gibbsianness is 
much more general and has nothing to do with microscopic reversibility. This can 
be checked readily for probabilistic cellular automata, [0, |^, but, more generally, it 
is the locality of the space-time interaction that does the job. In the present paper, 
we will exploit this fact in going from a global to a local fluctuation theorem for the 
entropy production in some models of interacting particle systems. At this moment, 
we need a second introduction to write about the statistical mechanics of steady state 
entropy production and how its fluctuations can give interesting information about 
the response of the system to perturbations. We refer to ^ for a recent impression. 
As we will see and as introduced in 0, once the conceptual framework and the Gibb- 
sian basis of the fluctuation theorem is understood, the transition from a global to 
a local fluctuation theorem will be merely a technical matter. Physically speaking 
however, it is much better for the obvious reason that global fluctuations are far too 
improbable to be observed, 0, [T^, |T^ . 
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1.1 Fluctuation theorem 



First observed in |]T3| and later derived in [O, IT5|, B^, the Gallavotti-Cohen fluctuation 



theorem proves a symmetry in the fluctuations in time of the phase space contraction 
rate for a class of dynamical systems. The dynamics must obey certain conditions; 
it is a reversible smooth dynamical system ^ i— >■ € i^'ona phase space K 

that is in some sense bounded carrying only a finite number of degrees of freedom (a 
compact and connected manifold). The transformation is a diffeomorphism on K. 
The resulting (discrete) time evolution is obtained by iteration and the reversibility 
means that there is a diffeomorphism 6 on K with 6'^ = 1 and 6 o (p o = (f)~^. It is 
assumed that the dynamical system satisfies some technical (ergodic) condition: it is 
a transitive Anosov system. This ensures that the system allows a Markov partition 
(and the representation via some symbolic dynamics) and the existence of the SRB 
measure p, an invariant measure with expectations 

1 ^ 



corresponding to time-averages for almost every randomly chosen initial point ^ & K 
(i.e., for an absolutely continuous measure with respect to the Riemann volume ele- 
ment on K) . The change of variables implied by the dynamics defines the Jacobian 
determinant J and one writes S = — In J. This is the phase space contraction rate 
which Gallavotti-Cohen identify with the entropy production rate via the following ar- 
gument: Define the (Shannon) entropy of a probability distribution m{d^) = m{^)d^ 
on K as 

S{m) = - j dim{i) lnm(0 (1.2) 
With m„ as density at time n, under the dynamics, the density at time + 1 is 

m^^M = (1.3) 
and the change in this entropy is therefore 



Af-l 



^ 

Dividing by and taking N to infinity, the empirical probability distribution ap- 
proaches the SRB distribution p, as in ( p..lD . Therefore, the time-averaged change 
in the entropy of the imagined reservoir is p{S), see also plj |23|. One further 
assumes (and sometimes proves) dissipativity: 

p{S) > 0. (1.5) 
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One is interested in the fluctuations of 



N/2 

in the state p. Informally, the fluctuation theorem then states that sn{0 has a 
distribution piy with respect to p such that 

lim^_ln^#, = l (1.7) 

^ Np{S)a Pn{-o) 

always. In other words, the distribution of the time-averaged S over long time in- 
tervals satisfles some general symmetry property. A more precise phrasing can be 
obtained via large deviation theory. For a continuous time version (Anosov flows) we 
refer to 0. 

The reason why we are interested in the fluctuation theorem is because the es- 
tablished symmetry in fluctuations is very general and it may be important for the 
construction of nonequilibrium statistical mechanics beyond linear order perturbation 
theory. Now there are various proposals, ideas and results but, at any rate, whatever 
the point of view, it is rather natural asking how to establish a local version of the 
fluctuation theorem. The title of the present paper refers to that problem with the 
understanding that local refers to a space-time window within a much larger spa- 
tially and temporally extended nonequilibrium system. This was already the subject 
of 1^, [TB[. We want to understand how general the local fluctuation theorem can be 
and what form it takes for some standard nonequilibrium models. It turns out that, 
as was explicitly pointed out in it is the Gibssian structure of the space-time 
distribution that allows a local fluctuation theorem. This was already apparent from 
1^, H but here we add further systematization and explain and illustrate this local 



version of the fluctuation theorem. 

One further question concerns the physical identiflcation of the quantity for which 
we are investigating the symmetry in the fluctuations. We will call it entropy pro- 
duction. This name already exists for a physical quantity that appears in close to 
equilibrium thermodynamics, and indeed we believe that our choice of words re- 
flects a generalization. The basic idea is that nonequilibrium steady states are not 
time-reversal invariant and that the mean entropy production should give a measure 
of discriminating between the original space-time distribution and its time-reversal. 
That is the relative space-time entropy density. For the variable entropy production, 
we must look up the source of the time-reversal symmetry breaking in the space-time 
interaction. It turns out that once it is recognized that the entropy production is the 
antisymmetric part of the space-time interaction under time-reversal, the symmetry 
in its local fluctuations (as expressed in the local fluctuation theorem (LFT)) is al- 
most an immediate consequence of the Gibbsian structure. This we will show. 
Of course, the question remains how we wish to use the local fluctuation theorem. 
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That is not the subject of the present paper but we refer to |T^, |T^, ^, ^, |T^, 
for some ideas. 



1.2 Example 

We sketch here the nature of a local versus global fluctuation theorem via a sim- 
ple model. We have in mind a (1 + l)-dimensional Ising spin system with formal 
Hamiltonian 

H{a) = J2M^)[(^t+i{x) + bat+i{x + l)] (1.8) 

x,t 

where we think x G Z as the spatial coordinate and t G Z as the (discrete) time; 
atix) = ±1, M 0. 
Look at the function 

T n 

Sn,t{(^) = ^ X] X] (^t{x)[at-i{x + 1) - crt+i{x + 1)] 

t=-T x=-n 

of the spins in a space-time window parametrized by ra, T > 0. We are interested in 
its fluctuations under the probability laws 

• Pn, the Gibbs measure on {—1, i}{-"'- -,n+i}x Z ^j^j-^ respect to the Hamiltonian 

n 

Hnicr) = 5Z crt{x)[at+i{x) + bat+i{x + 1)] 

t x=—n 

and 

■772 

• P, any infinite volume Gibbs measure on {—1,1}^ for the Hamiltonian ( |1.8| ). 

In both cases we take the counting measure as reference and set the inverse temper- 
ature P = 1. 

The difference is that P„ is an Ising model on a one- dimensional strip (finite spatial 
volume with infinite time-extension) and P is the corresponding model for infinite 
space-time volume. 

We start with the statement of a global fluctuation theorem; that concerns the law 
Pn- Consider the involution Qn,T by which all spins inside the window An^r = 
{—n, . . . , n + 1} X {— T — 1, . . . , T + 1} are reflected over the t = axis: (On,rC")i(a^) = 
<7-t{x) if {x,t) G An,T and remains unchanged otherwise. Remark that Sn,t{,Qn,TO') = 
—Sn,t{o') and upon writing Hn{Qn,TO') — Hn{cr) = 5'„,t(o") — i?n,T(c) we find, after a 
simple calculation, that for every function g 

I dPMgiSnA^)) = I rfP„(a)(7(-5„,T(^))e-^--('^)+^--('^) (1.9) 
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with \Bn^T{c)\ ^ en. 

As a result, for fixed n, for all functions /, 

■■■nil. I J''^-(--)/(g-.H--)/r) (1,0) 

which implies the symmetry expressed in (|1.7|) with N = T. 

Now to a local fluctuation theorem; that concerns the law P. A similar calculation 
shows that H{Qn,TO') — H{a) = Sn,t{o') — -B„,t((t) — Fn^ricr) with 

so that 

and we conclude that in both order of limits, 

n,TnT ^ j dP{a)f{-Sn,T{^)/W)e-^-'^^^y ^ ' 



This is the same symmetry as in (|1.71 ) but for the local fluctuations in a spatially 
extended system. Of course, ( |1.11| ) involves limits but the basic fact behind ( |1.11| ) is 



that there is a local function Rn,T = H o Qn,T — H, antisymmetric under the time- 
reversal Qn,T that preserves the a priori refeieiace measure, with \Rn,T{o') — Sn,T{,o')\ < 
Cin + C2T for which 



J dP{a)g{Rn,T) = J dP{a)g{-Rn,T)e- 



which is an exact local fluctuation symmetry. Various things are lacking from this 
example. Mathematically, things will be more complicated when the -Bn,T or Fn^r 
are not uniformly bounded or when time is not discrete or when the space-time 
Hamiltonian ( |1.8|) is not local or contains hard-core interactions. Physically, the 
example above carries no interpretation of Sn,T as entropy production. 



1.3 Local fluctuation theorem 

The main theme of the present paper is a general strategy to find a local fluctuation 
theorem for the entropy production in a nonequilibrium steady state, in the context 
of stochastic interacting particle systems. To get the idea we present the result infor- 
mally for a typical application. The details and mathematically precise statements 
about this model are given in Section 4. The model is a microscopic version of a 
react ion- diffusion system where the reaction consists of the birth and death of parti- 
cles on the sites of a regular lattice and the diffusion part lets these particles hop to 



6 



nearest neighbor vacancies subject to an external field. 

Consider the square lattice iB to each site i of which we assign a variable ri{i) = 0, 1, 
meaning that site is empty or occupied by a particle. The configuration ri can change 
in two ways: first, a particle can be created or destroyed at lattice site i: r\ ^ rf 
where rf is identical to 77 except that the occupation at the site i is flipped. Secondly, 
a particle at i can hop to one of the four nearest neighbor sites j under the condition 
that J is empty: 77 — ^ ri^^ where 77*-^ is the new configuration obtained by exchanging 
the occupations at sites i and j. We make a nonequilibrium dynamics by adding an 
external field E > ^ which introduces a bias for particle hopping in a certain direc- 
tion. 

In formula, first, a particle is destroyed or created at any given site at fixed rates. 
The transition from a configuration 77 to the new 77* takes place at rate 

c(«, n) = 7+(l - 77(i)) + 7_77(i) 

where 7+ is the rate for the transition — > 1 and 7_ is the rate for 1^0. Secondly, 
the particles on the lattice undergo a diffusive motion. To be specific, we choose a 
large square V centered around the origin with periodic boundary conditions and 
we first introduce hopping rates over a nearest neighbor pair (ij) in the horizontal 
direction, i = {11,12)-, j = (ii + 1,^2): 

c{i,j, ri) - e^/2^(z)(l - 77(j)) + e-^/^ri{j){l - ri{i)) 

The hopping rate in the vertical direction is constant (put £^ = in the above if 
j = (7i,72 i 1)). Taking E large, we expect to see many more jumps of particles 
to the right than to the left. In the absence of reaction rates, that is for 7± = 0, 
we recover the so called asymmetric exclusion process and particle number is strictly 
conserved. More generally, the Master Equation is 

= J2^PM)c{i,v') - Pt{v)c{i,v)] + J2^PM^)c{i,J,v'^) - Pt{v)c{i,j,v)] 

« {ij) 

For this model, the stationary measure p is the product measure with uniform density 
equal to 7+/ (7. + 7+) corresponding to a chemical potential ln7+/7_ of the particle 
reservoir. 

For a fixed nearest neighbor pair (ij), with j = {ii + 1,72) to the right of 7, the 
time-integrated microscopic current over an interval [— T, T] is 

with the total number of particles that have passed from site i to site j. We have 
the convention to take this current positive when the net number of particles jumping 
to the right (i.e., in the direction of the external field) is positive. Multiplying the 
sum of all the current contributions in V with the field E we get 

WvAVs, s e [-T, T]) = eY, 4{ris, s e [-T, T]) 

iev 
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which is a random variable representing the work done on our system over the time- 
interval [— T, T]. Its expectation in the stationary state equals (up to a temperature 
factor) the expected heat dissipated in the environment and is given by 

{Wv,t) = '2T\V\ Esmh(E/2) . ^+^- .„ 

(7+ + 7-J 

If we now fix another square A C inside our large system, then 

ieA:(n+l,i2)eA 

is the random variable "work done on the system in A over the time-interval [— T, T]" . 
That constitutes the main contribution to the local random variable "entropy produc- 
tion in the space-time window A x [— T, T]". Yet, this is only its bulk contribution. 
We have indeed only included in PVa^t the microscopic currents between the sites 
strictly inside A while particles will of course also hop in and out of A via its bound- 
ary. In other words, the region V\h. acts as a particle reservoir from which particles 
can enter or leave A. That also contributes to the entropy production as, quite gen- 
erally, the change in entropy in the particle reservoir equals the number of particles 
transferred to it, multiplied by its chemical potential. Now usually, this chemical po- 
tential is fixed and constant, i.e., not depending on whatever happens in the system 
itself. Here this is not the case. It suffices to imagine that almost all particles are 
in fact inside A with therefore a low density of particles in K \ A. As a result, the 
effective chemical potential for creating or destroying particles at the boundaries of A 
will depend on time and on whatever happened inside A before that time. Moreover 
this will contribute to the nonequilibrium condition only for E ^ because only then 
will there be a different rate of leaving/entering A at the right versus the left vertical 
boundaries of A. This is not the case for the upper versus the lower boundaries but 
also there, even when there would not be a field strictly inside A, the dynamics inside 
will be influenced by the flcld outside. This is summarized in the form of the second 
contribution to the time-integrated entropy production and it is a boundary term: 

<JdA,T = Re + Rr + Ru + Rd 

where the various terms correspond to the reactions taking place at the left, right, 
upper and lower boundaries of the square A. We will not write all of them down 
explicitly but here is for example 

where j = (ii -|- 1, 12), the sum over times t is over the times when a particle is created 
or destroyed at i, and the rates p and q are given by 

gA,,(j, ri, E) = e^/2prob[r^i(j) = 0|%(/c), keA,se [-T, t]] 
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and 



PA,t{j,V,E) = e-^/2Prob[r/,(j) = l\Vs{k),ke A, s G [-T,t]] 

where the probabihties refer to the steady state in V. In other words, the external 
field does not only work on the particles in A it also creates a gradient in chemical 
potential (large at the left boundary and smaller at the right) in A. The total random 
variable "entropy production in A" now reads 

Sa,T = Wa,T + J'dA,T 

The result proved in Section 4 is the fiuctuation theorem symmetry for Sa^t'- 



I Prob 



■Sa, 



T 



lim lim yyTTt; 111 — = a (1-12) 

."iAjr ^ 



A,T V \A\T Prob[^ 



uniformly in the j±. 

One may wonder whether the work W^a.t satisfies a similar fiuctuation symmetry. 
That is ( [ri2D with Wa,t repl acing Sa t- It remains uncertain however whether that 
is true uniformly in the values 7± i but, as we will show, it remains true whenever 

The rest of our paper is organized as follows: in Section 2 we give a general 
strategy to obtain LFT, which we apply in Section 3 for spinfiip processes and in 
Section 4 for the asymmetric exclusion process. 



2 Abstract setting 

We identify the essential mathematical structure, needed to pass from a global to a 
local fiuctuation theorem. Our later specific illustrations will then just be applica- 
tions of the same theme. 

We consider a measurable space {Q, JF) on which two sequences of probability mea- 
sures Pn and P^. Suppose that 9„ is an involution on Q such that P„ and Pn o 6„ 
are mutually absolutely continuous and the same for the pairs P^ and o 6„. We 
write 

/^^ = In , F„ = P„ + In 

dPA dP:: 

then, by definition, for all functions /, 

J dPJ{Rn) = j dP^e-^'-fi-Rr,) (2.13) 
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and 

J dPJ{R^) = J rfP„V-^"+^"/(-/?„) (2.14) 

The identity ( |2.13| ) expresses an exact symmetry in the fluctuations of i?„ but should 
be compared with the global symmetry (|1.9| , |1.10| ). The next equality (|2.14| ) is very 



similar but there is the correction term F„. To get rid of it (at least asymptotically 
in n) we need extra assumptions. This will then yield the local fluctuation theorem. 
Before we give a general way of expressing these assumptions, the reader may appre- 
ciate some more explication concerning our choice of 'global' versus 'local' as there 
is of course no natural interpretation of this within the proposed abstraction. 
As we will see in the next sections, we really start from two measures P and Pn on 
fl where Q will be the pathspace of an (infinite volume) interacting particle process 
on the (i— dimensional regular lattice Z*^; P will be an infinite volume steady state 
measure (i.e., the path-space measure of a stationary process over some time interval 
[— T, T]); n will refer to a finite space-time volume (corresponding to a sequence of 
cubes An centered around the origin times the interval [— T, T]) and G„ will be time- 
reversal on the space-time volume A„ x [— T, T]. The process P„ will be the path-space 
measure of the stationary interacting particle process on this finite A„ x [— T, T] . P^ 
is the marginal distribution of the trajectories restricted to the space-time window 
A„ X [— T, T] under P. In the context of interacting particle systems, P and P„ will 
be path-space measures of a Markovian process, whereas P^ will be non-Markovian. 
In the local fluctuation theorem it is attempted to recover the global symmetry of P„ 
under P„ also in the restrictions P^ of P to finite volumes A„. Clearly then, what we 
need is that the difference between P„ and P^ is a boundary term but this is more or 
less implied by having our interacting particle systems enjoy Gibbsianness on space- 
time. Finally, the meaning of P„ is that it gives, at least up to space-time boundary 
terms, a statistical mechanical representation of the thermodynamic steady-state en- 
tropy production. We wish however to refer to 0, ^ for explaining this. Still, 
it should be kept in mind that the P„ introduced in the following proposition will 
measure the difference between the true entropy production (denoted there by Sn) 
and Rn- 

There are in fact various strategies; we present two of them. 

Proposition 2.1: Let P„ be a measurable function so that B„ o 0„ = —B„. Define 
Sn = Rn + Bn and let (a„) be a sequence of positive numbers tending to infinity with 
n. Assume that P„ and P^ are mutually absolutely continuous and so that 

for all Ai, A2 G IR. Suppose that for all ^ G IR 

p(z)=lim — In / e-^^"dPn (2.16) 



n 
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exist and is finite. Then, whenever 



q{z) = hm — In / e-'^"dP'^ (2.17) 



n a 



n 



exists, then p{z) = q{z) and q{z) = q{l — z). 
Remarks: 

1. The symmetry q{z) = q{l — z) is dual to the symmetry as expressed in (|1.7|). Its 
Legendre transform i{a) = sup^{—q{z)—za) satisfies i{a) —i{—a) = —a. If Sn satisfies 
a large deviation principle under P„, respectively P^, then i{a) is the corresponding 
rate function, and the symmetry q{z) = q{l — z) is equivalent with the large deviation 
symmetry i{a) — i{—a) = —a. 

2. We will apply the strategy of Proposition 2.1 for obtaining a local fluctuation 
theorem for spinfiip processes in the next Section. 

Proof of Proposition 2.1: Since Bn o 0„ = —Bn, in the same way as for ( p.l3| ), 
we deduce that 

J dPJiSn) = J rfP„e-^"+^"/(-^„) (2.18) 

Starting with the left hand side, for /(s) = e"^**, by the Holder inequality, for 1/a + 
1/6=1 = l/v + l/w, 

In J dPne-'^" -l^""! "^^"^^^^ ^l^""! 

^lu/dP;(^)"-' (2.19) 



Dividing this by a„ and taking limits, we can use condition ( |2.15| ) with A2 = to get 



bv 

Again by the Holder inequality, both functions p and q are convex, and hence con- 
tinuous. Therefore we can take the limit for 6, f ^ 1 to conclude that p{z) = q{z). 



The right hand side of ( |2.18| ) can be treated in the same way: 



In [ rfP„e-(i-^)^"+^" <-\n[ ciP^(^)«e"^" + - In /" rfPV-''(^-^)^" 
J a J dP^ b J 

<-\n[ rfp;(^)'^e'^^" + I dPne-'"'^'-'^^-+''" + 

^In / rfP;(^)"'-ie-"'^"/" (2.20) 
bw J dPn 

11 



which, again after taking hmits n f +00, and using — — Rm gives 

qihil-z)) qi-hv{l-z) + l) 
p{z) = q{z) < ^ >> < ^ "--^ ^ 

and we can take the hmits b,v ^ 1 to get the desired q{z) — q{l — z) . 



Proposition 2.2: Let be a measurable function such that 5„ o ©„ = —Bn and 
define — Rn + B^- Let (a„) be a sequence of positive numbers tending to infinity 
with n so that for aU A e IR 

hm — In / dP^e^(-^"+^") = (2.21) 

Suppose that for aU 2; e IR 

q{z) = hm — In / e'^^'-dP^ (2.22) 

exists and is finite. Then, q{z) = q{l — z). 

Proof of Proposition 2.2: By definition of F„, we have 



We thus leave the left hand side and apply a similar chain of inequalities to the right 
hand side as was used in the proof of Proposition 2.1: 



-z)S„ 



< - In /" (iPre'^(^"+^") + — In /" dP'e-'''"^^-'^^-+^^+^^ + 
~ a J bv J 

ciP^e-"'^-^"-^")/'' (2.23) 

We may thus again divide by a„ and take limits first n t +00 to reach 

q^,) < - ^)) < g(-fa(l - ^) + 1) 

~ b ~ bv 
By convexity we can take the limits 6, J, 1 to obtain the desired conclusion. ■ 

Remarks: 

1. Of course, if it happens that |F„ + i?„|/a„ — > uniformly, then, for all positive 
functions /, 

hmlln^l™) ^0 
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without further ado. 

2. The difference between Proposition 2.1 and Proposition 2.2 is that in the first 
we suppose that Pn and are mutually absolutely continuous while in the latter, 
we need that P^ and -P^On are mutually absolutely continuous. We will follow the 
second strategy in Section 4 for the asymmetric exclusion process. 

3. The condition that the limits defining p{z) and q{z) exist is natural in the context 
where we have a large deviation principle for Sn under P^ and P^ resp. However if 
we define p'^ , p~ , , q~ by the corresponding limsup, resp. liminf, then we still have 
convexity of p"*", g"*" (the limsups), but not necessarily of p~ ,q~ ■ We can still conclude 
however the equality = q'^{z), and q~^{z) = q~^{l — z). 

3 LFT for spinflip processes 

We start our study with the, for physical applications, less interesting case of pure 
spinflip processes. For details on the construction of spinflip processes, we refer to 



The configuration space is-R' = {+l,— 1} (spins on the ci— dimensional regular lat- 
tice) and the path space is f2 = D{K, [— T, T]) the set of right-continuous trajectories 
having left limits, parametrized by time t G [—T, T], T > and having values ut G K. 
Our processes are specified in terms of spinflip rates c{x, a), x G Z'^, a G K for which 
our first most important assumption is that they are positive and bounded: there 
are constants bi > 0, 62 < +00 so that 61 < c{x, a) < 62 for all x, a. For convenience 
we assume that c{x, a) only depends on the neighboring spins (T{y) with ||/ — x| < 1. 
Thirdly, we assume the rates to be translation invariant: c{x,a) = c(0,r^(j). Here 
and afterwards we put A„ = [—n,nY H Z'^ 6„ denotes time-reversal on A„ defined 
by (0„a;)t(x) = uj-t{x) if x G A„, and (9„a;)t(a;) = ujt{x) if x ^ A„. On the jump- 
times we adapt Gnu; so that it becomes right- continuous, and thus obtain 9„ as an 
involution on Q 

We define A* = {x G A„, c{x, a) = c{x, a') for all a,a' E K with cr{y) = cr'{y), y G 
A„} for the subset of sites where the spinflip rates do not depend on the configuration 
outside A„. 

We first describe the sequence of processes P„ corresponding to P„ in the previous 
abstract setting. For this we fix a boundary condition r] & K and we define spinflip 
rates 

c„(x, a) = I[x G A„]c(x, aA„r/Ac ), X El^.a^K (3.24) 

where /[■] is the indicator function and (Jh^rj^c^ G K coincides with a on A„ and equals 
7] on the complement A^ = Z'^ \ A„. P„ is the stationary process on f2 with generator 

L„/(a) = 5^c„(x,a)[/(a^)-/(a)] 
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corresponding to a spinflip process in A„ with rates c{x, a) and boundary condition 
rj. We call the (unique) stationary measure p„ : / dpnL^f = 0. We always assume 
that for all a G {— 1, +1}^", p„((t) > 6i exp[— 62|A„|]. We can compute the density of 
Pn with respect to PnQn via a Girsanov formula for point processes, e.g. in p6| . 
For given uj & Q we let N^{uj), s G [—T, T], x G "E^ denote the number of spinflips at 
X up to time s; that is, N^{uj) = \{t G [—T,s],uJt-{x) = —ujt{x)}\] then, 

dPnBn J-T Cn{x,Us) Pn{uJT) 

As a consequence, the distribution of i?„ as induced from P„ satisfies immediately 
the global fluctuation symmetry ( |2.13| ). 
We can also consider the observable 



rp 



c{x,LJs-) 
c{x,UJs) 



which is measurable inside A„. This is not an arbitrary choice but there is too little 
physics here to call Rn or Sn the entropy production; we will not elaborate on this. 
The difference Bn = Sn — Rn is mainly a sum over a; G A„ \ A*. 
The other process P^ we need to look at is very similar but it is in general not Marko- 
vian. To define it, we take a stationary process P on Q and we take its restriction to 
An- We write p for the corresponding stationary measure and we let p„ be its restric- 
tion to An- We always assume that for all a G {—1, +1}^", Pn(o") > ai exp[— a2|A„|]. 
Being more explicit, we let P be an infinite volume stationary process with formal 
generator 

L/(a) = ^c(a;,a)[/K)-/(a)] 

X 

and put P^ the unique path-space measure such that 

1. The distribution of {ujt{x) : x G A„,t G [— T, T]} under P^ and P coincide. 

2. Under P^, ujt{x) = r]{x) for all x ^ A„, t>0. 
Theorem 3.1 [LFT for spinflip processes] For all z G IR, 

1 , f dPe~'^" 
lim —r- In = 

Proof of Theorem 3.1: 

Even though P^ is not Markovian (in general), it remains a jump process and the 
jump-intensities can be computed from the original spinflip rates. In order to have a 
Gibbsian structure these intensities must be the same in the bulk of A„ as they were 
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for the infinite volume process P. As the rates are local, the process restricted 
to An indeed has the same intensities as the process P„ except at the sites of the 
boundary A„, \ A*. This is a consequence of the following generally stated 
Lemma 3.2: 

Suppose Nt is a point process with intensity c^, i.e., Mt = Nt — Jq Cg ds is a martingale 
for the filtration JF^. Suppose that JF^' C J-* is a subfiltration of JF^, and define 



Then N'^ is a point process with intensity 



(3.25) 



(3.26) 



Proof of Lemma 3.2: 

It is easy to see that M[ = A*"^' — 1E(JJ Csds\J-'l.) = JEj[Mt\J^[\ is a J^j. martingale. Hence, 
it suffices to show that 



Bt = JE[[ Csds\J='l] - [ W:[cs\J^'^]ds 
Jo Jo 



is a jF^'-martingale. This is a consequence of the following equalities: 

= R + 1E 



Cr - ^[Cr\J^',])dr\r^ 



(3.27) 



(3.28) 



Therefore, the rates of the restricted process on A„ are given by 

Cs{x,uj) = lE[c(x,o-,)|cr^(y) = uJriy), -T <T < s,y E A„] 
where the expectation is with respect to P. 

Or, for all x G A„, Nf{uj) — J^Cs{x,u) is a martingale under P^. As a consequence, 
Cs{x, u) = c{x, a) when x G A* and uig = cr. 

Just as for the pair P„, PnQn, the absolutely continuity of P^On with respect to P^ and 
vice versa is guaranteed by the positivity of the spinflip rates inside A„. We are thus 
ready to apply Proposition 2.1. We must first verify the corresponding assumption 
(glSl) . We find 



Bn Sri. Rn 



J2 jdNy{u)\n 



a- 6 An \ A 



(3.29) 
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and 



dPn Pn[0J-T) ^/r^K^J Cn[x,UJs) J -T 

(3.30) 



xgA„\A* 



Clearly, both \Bn\ and | ln(iP^/(iP„| are bounded by ciA^([-T, T], A„ \ A^) + C2|A„| + 
C3T|A„ \ A* I for some constants ci,C2,C3 < oo, where iV([— T, T], A„, \ A*) is the 
number of spinflips that have occurred in the space-time window [— T, T] x (A„ \ A*). 
It remains thus to show for all A 

lim^ln [ rfp^V^([-r-riAAA*) 

nToo n° J 

= hm ^ In /" rfp0^gA7V([-T,T],A„\A;) ^ 

where we inserted the reference process P*^ (and its restriction P° to A„) correspond- 
ing to the product process of independent spinflips (rate 1). In particular, 

hm A In / rfpOe^^d-^'^l'^"^'^") = hm ^^A^2Tre^ -11=0 



ntoo n'^ J ntoo n'^ 



and we can apply the same argument as in Proposition 2.1. Finally, the condition 
(p.l6|) of Proposition 2.1 is a consequence of the large deviation results of [p4i. ■ 



4 LFT for the asymmetric exclusion process 

The configuration space is now K = {0, 1} (occupation variables on the 2-dimensional 
regular lattice) and the pathspace Q = D{K, [— T, T]) is essentially unchanged from 
that in the previous section. For t] G K, ri{x) = 1,0 indicates the presence, respec- 
tively absence of a particle at the site x . This hopping dynamics will be modeled 
by an asymmetric exclusion process. This is a bulk driven diffusive lattice gas. The 
hopping rates for vertical (f) and horizontal bonds iK) depend on the direction in 
the following way: 

c''(x,?7) = \Aj]{x){\ - r]{x \ e-i)) +?7(x + e2)(l -?7(a;))] 

c'^'^(x, ri) = ^[e^/2^(x)(l - 7]{x + d)) + e~^/^r]{x + ei)(l - r]{x))] 

where ei, 62 are the unit vectors in the positive horizontal and vertical direction. 
In addition, for the moment, we allow for the possibility of particle creation and 
destruction. We put the birth/death rate c{x, rf) = e independent of the configuration 
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rj and the site x. 

The formal Markov generator L to the infinite volume process is then found as the 
sum 

X (xy) 

where 77^ is the new configuration after changing the occupation at x, t]^"^ is the 
new configuration after exchanging the occupations at x and y and c(x, y, rj) is given 
by ( [4.31|) for nearest neighbors x,y = x ± ei. We can allow more general reaction- 
diffusion processes (e.g. with extra interaction, speed change, etc.) but we will stick 
here to this choice. What is simpler here is that the Bernoulli measure p with density 
1/2 is a non-reversible stationary measure. The corresponding pathspace measure 
over the time-interval [— T, T] is P = and we put the process restricted to 
the finite square A„ This P^ will now play the role of P^ of Section 2. For a given 
trajectory uj E Q we let N^y{uj), s G [— T, T] be the number of hopping times where 
the occupation at the nearest neighbor sites (xy) was exchanged. Since this model 
has a clear physical interpretation we can define the variable entropy production in 

An- 

The first contribution comes from the work done by the external field 

Wn = E Yl f dN:y{u)[us-{x){l-uJs~{y))-uj,^{y){l-uj,-{x))] 

(4.31) 



(xj/),j;=x+eiGA„ ' 



This is of the form field {E) times current. There is a second contribution from 
differences in the reaction rates at the boundary: particles enter or leave at different 
rates at the various boundaries; this contribution is present even in the case where 
no external field is applied inside A^: 



2e + A^(^,s)e^^/2 
^ . j-T 2e + A-^(cj,s)e-^^/2 



2e + «:^(a;,s)e-^^/2 



+ {l-ujs-{x))\n- 



2e + K-^(cj,s)e^j/2^ 



2e + K^{uj,s) 



where the second sum is over all (external) neighbors y of x and dAn is the interior 
boundary of A„. Here, the additional rates are 

Af (r/, t) = 1E^[1 - 7]tiy) \r],iz), z e An, s e [-T, t]] 
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and 

K^ir],t) = JE''[r]t{y)\risiz),ze A„, s G [-T,t]] 

for Ej = ±E ii y = X ± ei, Ej = if y = X ± 62 and the expectations are in the 
process P = P^. The variable entropy production is put 

Sn = Wn + Jn 

The symmetry in the fluctuations of Sn is given by 

Theorem 4.1 [LFT for the aymmetric exclusion process] For all e (including e = 0), 
for all z ^ R, 

1 , f dPe-'^" 
lim In — ^ = 

n,T n^T J dPe-^^-^)^- 

Proof of Theorem 4.1 We start by noting that for the time-reversal 0„, 

pEo ^ p-E 

Obviously then, for a function / measurable in A„ x [— T, T], 

dPfiu:) = I dP^fico) = I dP^^m^u) 

and we must investigate express the density dP^ / dP~^ via a Girsanov formula. This 
is the strategy of Proposition 2.2. The Girsanov formula gives 

dP^ 

^^^ = Sn + F^ (4.33) 
with the following correction term: 

F„H = sinh(E/2) / ds[u;s{y){l-u;s{x))-u;s{x){l-u;,{y))] 



+ E E / ds[u;s{x)[X^{uJ,s))-X-'^{u;,s)] 

xGdA y=x:Lei ,xdie2GA^ 

+ (1 - uJs{x))[k^{uj, s)) - s)] (4.34) 



Now, < c\dA\T because the (first) bulk term in ( |4.34|) telescopes to a boundary 
term. We can thus apply Remark 1 after Proposition 2.2 to finish the proof. ■ 

Next, we investigate whether the variable work Wn of ( [4.3 1| ) itself satisfies the 
same local symmetry as the entropy production. 
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Theorem 4.2 [LFT for the work done] For ah e > 0, for aU z G H, 

1 JdPe-'^" 

Proof of Theorem 4.2 Clearly, since e > 0, \Jn\ < cA^([-r, T], A„ \ A;), that is 
bounded, up to a constant, by the number of flips in the trajectory on sites x G A„\A* , 
for times t G [— T, T]. We can therefore verify condition ( p.21| ) in the same way as 
we did for Theorem 3.1. Finally, the large deviation results of remain valid for 
e > 0, so that we can finish the proof along the lines of Proposition 2.2. ■ 



5 Remarks 

1. It is clear from the preceding analysis that the reasons for having a global 
or local fluctuation theorem do not in any way depend on the 0„ being time- 
reversal. Thus, the same results will be reproduced in exactly the same form for 
any other involution. Of course, the symmetry breaking part in the pathspace 
action functional will be the variable for which the fluctuation symmetry holds 
(replacing entropy production corresponding to time-reversal symmetry break- 
ing). As an example, if for a spinflip process, the rates are not even under a 
global spinflip (by the presence of a bias or magnetic field), then a local fluctu- 
ation theorem will be established for the variable magnetization. Furthermore, 
we may consider the composition of two or more involutions — in this way, we 
could e.g. obtain a local fluctuation theorem for the odd part (under spinflip) 
of the variable entropy production. Finally, we can even go beyond the case 
of involutions and consider instead the generators of the symmetry group for 
the unperturbed dynamics. In this case, the precise form of the fluctuation 
symmetry is not preserved but its modification presents no real problem. 

2. We restricted our discussion to interacting particle systems where the evolution 
is Markovian. Within the Gibbsian space-time picture, this means that the 
interaction is "nearest neighbor" in the time direction (the jump intensity at 
time t depends only on the configuration at time t~). However, this restriction is 
not at all necessary. If the jump intensities are local in space and bounded from 
above and from below, then we can still apply the Girsanov formula for point 
processes to obtain the local fluctuation theorem from the global fluctuation 
theorem. 
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